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C\j ! Abstract 
> 

We generalize the quasilocal definition of the stress energy tensor of Einstein gravity to the 
case of third order Lovelock gravity, by introducing the surface terms that make the action well- 



defined. We also introduce the boundary counterterm that removes the divergences of the action 
and the conserved quantities of the solutions of third order Lovelock gravity with zero curvature 
boundary at constant t and r. Then, we compute the charged rotating solutions of this theory 
in n + 1 dimensions with a complete set of allowed rotation parameters. These charged rotating 
solutions present black hole solutions with two inner and outer event horizons, extreme black 
holes or naked singularities provided the parameters of the solutions are chosen suitable. We 
compute temperature, entropy, charge, electric potential, mass and angular momenta of the black 
hole solutions, and find that these quantities satisfy the first law of thermodynamics. We find 
a Smarr-type formula and perform a stability analysis by computing the heat capacity and the 
determinant of Hessian matrix of mass with respect to its thermodynamic variables in both the 
canonical and the grand-canonical ensembles, and show that the system is thermally stable. This 
is commensurate with the fact that there is no Hawking-Page phase transition for black objects 
with zero curvature horizon. 
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I. INTRODUCTION 



In four dimensions, the Einstein tensor is the only conserved symmetric tensor that de- 
pends on the metric and its derivatives up to second order. However for spacetimes possessing 
more than four dimensions, as assumed in both string theory and brane world cosmology, 
this is not the case. In string theory, extra dimensions are a theoretical necessity since 
superstring theory requires a ten-dimensional spacetime to be consistent from the quantum 
point of view, while in brane world cosmology matter and gauge interactions are localized 
on a 3-brane, embedded into a higher dimensional spacetime in which gravity propagates 
throughout the whole of spacetime. The most natural extension of general relativity in 
higher dimensional spacetimes with the assumption of Einstein - that the left hand side 
of the field equations is the most general symmetric conserved tensor containing no more 
than two derivatives of the metric - is Lovelock theory. Lovelock [1] found the most general 
symmetric conserved tensor satisfying this property. The resultant tensor is nonlinear in 
the Riemann tensor and differs from the Einstein tensor only if the spacetime has more 
than 4 dimensions. Since the Lovelock tensor contains metric derivatives no higher second 
order, the quantization of the linearized Lovelock theory is ghost-free The concepts of 
action and energy-momentum play central roles in gravity. However there is no good local 
notion of energy for a gravitating system. Quasilocal definitions of energy and conserved 
quantities for Einstein gravity 3, 0| define a stress energy tensor on the boundary of 
some region within the spacetime through the use of the well-defined gravitational action 



of Einstein gravity with the surface tern) of Gibbons and Hawking (J Our first aim in 
this paper is to generalize the definition of the quasilocal stress energy tensor for computing 
the conserved quantities of a solution of third order Lovelock gravity with zero curvature 
boundary. The first step is to find the surface terms for the action of third order Lovelock 
gravity that make the action well-defined. These surface terms were introduced by Myers 
in terms of differential forms The explicit form of the surface terms for second order 
Lovelock gravity has been written in Ref. ja|. Here, we write down the tensorial form of the 
surface term for the third order Lovelock gravity, and then obtain the stress energy tensor 
via the quasilocal formalism. Of course, as in the case of Einstein gravity, the action and 
conserved quantities diverge when the boundary goes to infinity. We will also introduce a 
counterterm to deal with these divergences. This is quite straightforward for the cases we 
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consider in which the boundary is flat. This is because all curvature invariants are zero 
except for a constant, and so the only possible boundary counterterm is one proportional 
to the volume of the boundary regardless of the number of dimensions. The coefficient of 
this volume counterterm is the same for solutions with flat or curved boundary. The issue 
of determination of boundary counterterms with their coefficients for higher-order Lovelock 
theories is at this point an open question. Since the Lovelock Lagrangian appears in the low 
energy limit of string theory, there has in recent years been a renewed interest in Lovelock 
gravity. In particular, exact static spherically symmetric black hole solutions of the Gauss- 
Bonnet gravity (quadratic in the Riemann tensor) have been found in Ref. [jj, and of the 
Maxwell-Gauss-Bonnet and Born-Infeld-Gauss-Bonnet models in Ref. ^(|. The thermody- 



namics of the uncharged static spherically black hole solutions has been considered in 



of solutions wit 
solutions in 10, 



i uoutrivial topology and asymptotically do Sitter in Q and of charged 



111, 



13J. Very recently NUT charged black 



no 



e solutions of Gauss-Bonnet grav- 



ity and Gauss-Bonnet-Maxwell gravity were obtained l^J. All of these known solutions in 
Gauss-Bonnet gravity are static. Not long ago one of us introduced two new classes of rotat- 
ing solutions of second order Lovelock gravity and investigated their thermodynamics jl^j ]. 
and made the first attempt for finding exact solutions in third order Lovelock gravity with 
the quartic terms ^(|. Our second aim in this paper is to obtain rotating asymptotically 
anti de Sitter (AdS) black holes of third order Lovelock gravity and investigate their ther- 
modynamics. Apart from their possible relevance to string theory, it is of general interest 
to explore black holes in generalized gravity theories in order to discover which properties 
are peculiar to Einstein's gravity, and which are robust features of all generally covariant 
theories of gravity. The outline of our paper is as follows. We give a brief review of the 
field equations of third order Lovelock gravity and the counterterm method for calculating 
conserved quantities in Sec. HTJ In Sec. IIHI we introduce the (n + l)-dimensional solutions 
with a complete set of rotational parameters and investigate their properties. In Sec. IIVI 
we obtain mass, angular momentum, entropy, temperature, charge, and electric potential 
of the (n + l)-dimensional black hole solutions and show that these quantities satisfy the 
first law of thermodynamics. We also perform a local stability analysis of the black holes 
in the canonical and grand canonical ensembles. We finish our paper with some concluding 
remarks. 
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II. FIELD EQUATIONS 



The action of third order Lovelock gravity in the presence of electromagnetic field may 
be written as 

Ig=^t I d n+1 x^g~ (-2A + R + a 2 C 2 + a 3 C 3 - F^F^) (1) 
lovr J M 

where A is the cosmological constant, a 2 and 0:3 are Gauss-Bonnet and third order Lovelock 
coefficients, F pv = d p A v — d v A p is electromagnetic tensor field and A p is the vector potential. 
The first term is the cosmological term, the second term, R, is the Einstein term, the third 
term is the Gauss-Bonnet Lagrangian given as 

C 2 = R^ S R^ 5 - ±R^R» V + R 2 (2) 
and the last term is the third order Lovelock term 

C O Tyfll/aK p DPT I Q TDfJ-f T3 aK T3P T I O/l DM^ CTK P PP tO\ 

L 3 — ArT ix (JK p T tx ^ u + o-ix up-ri VT ix M(t + l^HT ti aRvp tir p [6 J 

+3RR >ll " 7K R aKpu + 2AR llV(IK R up R KV + IQR^ R va R a M — 12RW 1V 'R^ + R 3 

From a geometric point of view the combination of these terms in seven and eight dimen- 
sions is the most general Lagrangian that yields second order field equations, as in the 
four- dimensional case for which the Einstein-Hilbert action is the most general Lagrangian 
producing second order field equations, or the five- and six-dimensional cases, for which the 
Einstein-Gauss-Bonnet Lagrangian is the most general one fulfilling this criterion. Since the 
third Lovelock term in eq. (0) is an Euler density in six dimensions and has no contribution 
to the field equations in six or less dimensional spacetimes, we therefore consider (n + 1)- 
dimensional spacetimes with n > 6. Varying the action with respect to the metric tensor 
g pu and electromagnetic tensor field F pv the equations of gravitation and electromagnetic 
fields are obtained as: 

Gjg + A 9lw + a 2 G$ + a 3 G$ = T„ v (4) 
VuF^ = (5) 

where T pu = 2F p ^F pu — \F pa F p(J g pv is the energy-momentum tensor of electromagnetic field, 
G$ is the Einstein tensor, and Gfi, and Gfi, are the second and third order Lovelock tensors 
given as jl7|: 

= 2{R paiiT RJ KT — 2R ppua R fXT — 2R pa R° v + RR pu ) — -C-zQuv (6) 
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/^(3) Q f/l P r P°" K P pA QP t P dtk pA I o d ran td pAp 

~~ -<H 4 -K ^ K aK \ p K VTil - oK Xa R Tfl K vpK + ZR U R aK \ p R TjU 

ri ' rL GKTp rL VfJ, T O.TL Uap -TL rp 171 K ' 0rt UTK 11 up* 1 p 
I A p TCTK p pp p T<TK p pp _|_ TJTpuK p p _|_ O P P KT P P 

+8R T uppR p a RCT T - %R a vT P R T oRp - 8R T apR c ' T Rv P - 4:RR T I/fip R P T 
+AR rp R pr R up — 8R T u R rp R p p + 4RR up R p fl — R 2 R up ) — -C^g^ (7) 

The Einstein-Hilbert action (with «2 = «3 = 0) does not have a well-defined variational 
principle, since one encounters a total derivative that produces a surface integral involving 
the derivative of 5g pu normal to the boundary. These normal derivative terms do not vanish 
by themselves, but are canceled by the variation of the Gibbons-Hawking surface term |f| 

4 X) = lzl d n x^K (8) 
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The main difference between higher derivative gravity and Einstein gravity is that the surface 
term that renders the variational principle well-behaved is much more complicated. However, 
the surface terms that make the variational principle well-defined are known for the case of 
Ga^BonnetgravityBatobe^'^f'whe^f'is 

and where 7^ is induced metric on the boundary, K is trace of extrinsic curvature of 
boundary, G^ b is the n-dimensional Einstein tensor of the metric 7 a b and J is the trace of 

J ab = l -{2KK ac K c b + K cd K cd K ab - 2K ac K cd K db - K 2 K ah ) (10) 



For the case of third order Lovelock gravity, the surface term that makes the variational 
principle well defined is I b = I b + I b + I b , where 1^ is 

4 3) = [ d n x^{3a 3 (P - 2G^K ab - 12R ab J ab + 2RJ 

JSM 

-AKR abcd K ac K M - 8R abcd K ac K b K ed )} (11) 

In eq. fTTf G a 2 b is the second order Lovelock tensor (JHJ) for the boundary metric ^ ab , and P 
is the trace of 

p ab = ^{[ir 4 -6if 2 ir d K cd + 8i^^ 

-(AK 3 - 12KK ed K ed + %K de K)K^ d )K ac Kl - 2AKK ac K cd K de K e b 

+ {12K 2 - 12K ef K £ f)K ac K cd K db + 2AK ac K cd K de K e f K bf } (12) 



In general Iq+I^ + is divergent when evaluated on solutions, as is the Hamiltonian 



and other associated conserved quantities 
is through the use of background subtraction 



. One way of eliminating these divergences 
in which the boundary surface is embedded 



in another (background) spacetime, and all quasilocal quantities are computed with respect 
to this background, incorporated into the theory by adding to the action the extrinsic cur- 
vature of the embedded surface. Such a procedure causes the resulting physical quantities 
to depend on the choice of reference background; furthermore, it is not possible in general 
to embed the boundary surface into a background spacetime. For asymptotically AdS so- 
lutions, one can instead deal with these divergences via the counterterm method inspired 
by AdS/CFT correspondence ^|. This conjecture, which relates the low energy limit of 
string theory in asymptotically anti de-Sitter spacetime and the quantum field theory on its 
boundary, has attracted a great deal of attention in recent years. The equivalence between 
the two formulations means that, at least in principle, one can obtain complete information 
on one side of the duality by performing computation on the other side. A dictionary trans- 
lating between different quantities in the bulk gravity theory and their counterparts on the 
boundary has emerged, including the partition functions of both theories. In the present 
context this conjecture furnishes a means for calculating the action and conserved quantities 



intrinsically without reliance on any reference spacetime |19l I20L |21| by adding additional 
terms on the boundary that are curvature invariants of the induced metric. Although there 
may exist a very large number of possible invariants one could add in a given dimension, 
only a finite number of them are nonvanishing as the boundary is taken to infinity. Its 
many applications include computations of conserved quantities for black holes with rota- 
tion, NUT charge, various topologies, rotating black strings with zero curvature horizons 
and rotating higher genus black branes 2^|. Although the counterterm method applies for 
the case of a specially infinite boundary, it was also employed for the computation of the 
conserved and thermodynamic quantities in the case of a finite boundary 23[. Extensions 
to de Sitter spacetime and asymptotically flat spacetimes have also been proposed Q- All 
of the work mentioned in the previous paragraph was limited to Einstein gravity. Here we 
apply the counterterm method to the case of the solutions of the field equations of third 
order Lovelock gravity. At any given dimension there are only finitely many counterterms 
that one can write down that do not vanish at infinity. This does not depend upon what 
the bulk theory is - i.e. whether or not it is Einstein, Gauss-Bonnet, 3rd order Lovelock, 
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etc. Indeed, for asymptotically (A)dS solutions, the boundary counterterms that cancel 
divergences in Einstein Gravity should also cancel divergences in 2nd and 3rd order Love- 
lock gravity. The coefficients will be different and depends on A and Lovelock coefficients 
as we will see this for the volume term in the flat boundary case below. Of course these 
coefficients should reduce to those in Einstein gravity as one may expect. Unfortunately we 
do not have a rotating solution to either Gauss-Bonnet or 3rd-order Lovelock gravity that 
does not have a flat boundary at infinity. Consequently we restrict our considerations to 
counterterms for the flat-boundary case, i.e. R a bcdi.l) — 0, for which there exists only one 
boundary counterterm 

I« = ^-t d^V^^ 1 , (13) 
8tt Jsm l 

where L is a scale length factor that depends on I, a 2 and 013, that must reduce to I as a 2 
and a 3 go to zero. Having the total finite action I = I G + 1^ + 1^ + 1^, one can use the 
quasilocal definition P, |J to construct a divergence free stress-energy tensor. For the case 
of manifolds with zero curvature boundary the finite stress energy tensor is 

T ab = —{{K ab - K-f ab ) + 2a 2 {3J ab - J~f ab ) 

nn 1 

+3a 3 {5P ab - PY b ) + —^l ah }• (14) 

The first three terms in eq. (fTlj) result from the variation of the surface action (|SJ)-([12)) with 
respect to 7 afe , and the last term is the counterterm that is the variation of I ct with respect 
to 7 ab . To compute the conserved charges of the spacetime, we choose a spacelike surface B 
in dM. with metric aij, and write the boundary metric in ADM form: 

l ab dx a dx a = -N 2 dt 2 + (Tij (dtp* + V l dt) (dtpP + V j dt) , (15) 

where the coordinates ip l are the angular variables parameterizing the hypersurface of con- 
stant r around the origin, and N and V 1 are the lapse and shift functions respectively. When 
there is a Killing vector field £ on the boundary, then the quasilocal conserved quantities 
associated with the stress tensors of eq. ((Hj) can be written as 

Q(0 = / d^V^T^e, (16) 

JB 

where a is the determinant of the metric a^, and n a is the timelike unit normal vector to 
the boundary B. For boundaries with timelike (£ = d/dt) and rotational (q = d/d<p) Killing 
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vector fields, one obtains the quasilocal mass and angular momentum 

M = f d n -\^T ah n a i\ (17) 
Jb 

J = [ d n '\^T ah n\\ (18) 
Jb 

provided the surface B contains the orbits of These quantities are, respectively, the 
conserved mass and angular momentum of the system enclosed by the boundary B. Note 
that they will both depend on the location of the boundary B in the spacetime, although 
each is independent of the particular choice of foliation B within the surface dJ\A. 



III. (n + l)-DIMENSIONAL ROTATING SOLUTIONS 

As stated before, the third order Lovelock term in eq. (0) is an Euler density in six 
dimensions and has no contribution to the field equations in spacetimes of dimension six 
or less. Taking n > 6, we obtain the (n + l)-dimensional solutions of third order Lovelock 
gravity with nonvanishing electromagnetic field with k rotation parameters and investigate 
their properties. The rotation group in n + 1 dimensions is SO(n) and therefore the number 
of independent rotation parameters is [(n + l)/2], where [x] is the integer part of x. The 
metric of an (n + l)-dimensional asymptotically AdS rotating solution with k < [(n + l)/2] 
rotation parameters whose constant (t. r) hypersurface has zero curvature may be written 
as 3 

/ k \ 2 r 2* _ , 2 

ds 2 = —f(r) | Edt — '^^a i d(p i J + -j^^ (a-idt — El 2 d(j)ij 



I 4 

t=i / i=i 



■— - - — 22(<nd<f>3 ~ a^i? + r 2 dX 2 , (19) 

J^ r ! i<j 



where H = y 1 + Yli a i/l 2 i ^e angular coordinates are in the range < <pi < 2n and dX 2 
is the Euclidean metric on the (n — k — l)-dimensional submanifold with volume E n _fc_i. 
Using eq. (|5)i. one can show that the vector potential can be written as 



A »= {n -l)rn-2 ^ 5 l- a i 5 % (no sum on i). (20) 



where q is an arbitrary real constant which is related to the charge of the solution. To find 
the function f(r), one may use any components of eq. (@J). The simplest equation is the rr 
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component of these equations which can be written as 

" ' /' + Ar 6 

+360(VW 3 - 12(V)aarV a + CV V/ = -q 2 r 8 ~ 2n (21) 



T) — 1 



where prime denotes the derivative with respect to r. We present the full solution to this 
equation in the appendix. Here, for simplicity, we consider the solutions of eq. (|21j) for a 
restricted version of the aj's given as 



a 2 <\ 



" 3 = 72C7 1 ) ' = (n-2)(n-3) (22) 

Equation (J2"T|) with condition (}2*2*|) has one real and two complex solutions that are the 
complex conjugate of each other. The real solution of eq. (J2*T|) with condition is 

1 ' 7 6Act 3am 6ag 2 \ 1 ^ 3 . 

+ ~ 17 71 7771 o\ „9r,.-9 r (23) 




n(n — 1) r n (n — 1) (n — 2)r 2n 2 

where m is the mass parameter. Although the other components of the field equation (0J) 
are more complicated, one can check that the metric (jl9j) satisfies all the components eq. (j3J) 
provided above f(r) is given by (|23jl . Unlike the solution in Gauss-Bonnet gravity, which has 
two branches, the solution (|2~3j) has only one branch. Indeed, eq. (f2*Tj) with the above a's has 
the real solution (|2~H|) and two complex solutions that are complex conjugates of each other. 
This feature is the same as the asymptotically AdS solution of Ref . |26( , which has a unique 
anti de Sitter vacuum. These solutions are asymptotically AdS or dS for negative or positive 
values of A respectively. In this paper we are interested in the case of asymptotically AdS 
solutions, and therefore we put A = — n(n — l)/2l 2 . One can show that the Kretschmann 
scalar R^ v \ k R^ uXk diverges at r = 0, and therefore there is a curvature singularity located at 
r = 0. Seeking possible black hole solutions, we turn to looking for the existence of horizons. 
As in the case of rotating black hole solutions of Einstein gravity, the above metric given by 
eqs. (|T^|) and ()46|) has both Killing and event horizons. The Killing horizon is a null surface 
whose null generators are tangent to a Killing field. The proof that a stationary black hole 
event horizon must be a Killing horizon in the four- dimensional Einstein gravity [22( cannot 
obviously be generalized to higher order gravity. However the result is true for all known 
static solutions. Although our solution is not static, the Killing vector 

k 

X = d t + YP^ ( 24 ) 



»=i 



is the null generator of the event horizon, where k denotes the number of rotation parameters. 
The event horizon is defined by the solution of g rr = f(r) = 0. For the case of uncharged 
solutions, there exists only one event horizon located at 

r+ = (m/ 2 ) 1/n (25) 

This feature is different from the case of unc harg ed spherically symmetric solutions of third 
order Lovelock gravity with curved horizon [lq, for which one may have uncharged black 
holes with two horizons, extreme ones, or a naked singularity. As we demonstrate in the 
appendix, the general uncharged solution does not have two horizons. The charged solution 
presents a black hole solution with two inner and outer horizons, provided the mass param- 
eter m is greater than, m ext , an extreme black hole for m = m ext and a naked singularity 
otherwise, where m ext is 

2(n-l) / 2q 2 l 2 X"/ 2 ^- 1 ) 

m ^ = ^2jp{^T)) (26) 

The general charged solution also has only two horizons, as shown in the appendix. 

IV. THERMODYNAMICS OF BLACK HOLES 

One can obtain the temperature and angular momentum of the event horizon by analytic 
continuation of the metric. Setting t — *■ ir and a% — > iai yields the Euclidean section of (|19jh 
whose regularity at r = r + requires that we should identify r ~ r + j3 + and fa ~ fa + 
where (3 + and Q^s are the inverse Hawking temperature and the angular velocities of the 
outer event horizon. One obtains 

o-x _ m _ f'(r + ) _ n(n-l)-2lYr- 2(n ~ 1) 

a = (28) 

Next, we calculate the electric charge of the solutions. To determine the electric field we 
should consider the projections of the electromagnetic field tensor on special hypersurfaces. 
The normal to such hypersurfaces is 

u° = ^, u r = 0y = ^, (29) 
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and the electric field is = g^ p F pv u v . Then the electric charge per unit volume V n -\ can 
be found by calculating the flux of the electric field at infinity, yielding 

Q=h q (30) 



The electric potential measured at infinity with respect to the horizon, is defined by 33 1 



$ = a^Uco - A^U + > (3i) 



where \ is the null generator of the horizon given by eq. (|24jl . We find 



* = 7 (32) 

Conserved quantities associated with the spacetime described by (fT9*j) can be obtained via 
the counterterm method. Using eqs. ()17|) and (jl8j) . the mass and angular momentum will 
be finite provided 

l5l 2 y/a{l-\) 



L 



hi 2 + 9a - l 2 \ 2 - Al 2 X 



Q \ 1/3 
601 x 



I 2 

where we note that L reduces to I as a goes to zero. The mass and angular momentum per 
unit volume V n -\ can then be obtained through the use of eqs. (fTTjl and (fT£j) . We find 

M = -^[nH 2 -l]m, (33) 

l07T 

Jj = nEmcii (34) 

16tt K J 

Black hole entropy typically satisfies the so-called area law, which states that the entropy 
of a black hole equals one-quarter of the area of its horizon This near-universal law 
applies to all kinds of black holes and black strings in Einste in g ravity j^j. However in 



higher derivative gravity the area law is not satisfied in general 



gravity that 



31 



32] 



It is known in Lovelock 



[(d-l)/2] 



S = \ E kak I d n - l x^gC k ^ (35) 
k=l ^ 

where the integration is done on the (n — l)-dimensional spacelike hypersurface of Killing 
horizon, g^ v is the induced metric on it, g is the determinant of g^ and Ck is the kth order 
Lovelock Lagrangian of g pv . For the topological class of black holes we are considering, the 
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horizon curvature is zero. Consequently the area law holds. Denoting the volume of the 
hypersurface boundary at constant t and r by V n -\ = (27r) fc £ n _jfc_i, we obtain 

S = =Vr 1 (36) 

for the entropy per unit volume V n -\. The entropy can also be obtained through the use of 
Gibbs-Duhem relation 

S = /3{M - T t C t ) - I (37) 

where I is the finite total action evaluated on the classical solution, and d and r« are the 
conserved charges and their associate chemical potentials respectively. For simplicity, we 
consider the uncharged solutions, for which Q = Ji and Tj = fij. Using eqs. ((H), (jEJ-JHJ) 
and (JT3J), the finite total action per unit volume V n ~\ can be calculated as 

I = —^rl (38) 

Now using Gibbs-Duhem relation (J37)l and eqs. andlB^j) one may confirm that the 

entropy per unit volume obeys the area law of eq. ([36)1 . 

A. Energy as a function of entropy, angular momenta and charge 

We first obtain the mass as a function of the extensive quantities S, J and Q. Using 
the expression for the entropy, the mass, the angular momenta, and the charge given in eqs. 
(|27jl . Q3Up. and (|34p. and the fact that f(r + ) = 0, one can obtain a Smarr-type formula 

as 

where J 2 = | J| 2 = Jf and Z = H 2 is the positive real root of the following equation 

(Z - l)^ - — { Mn-l)(n-2)lSJ 1^ 

1 ; 165 2 \n[(n-l)(n-2)5 2 + 27r 2 g 2 / 2 ] J 1 ; 

One may then regard the parameters S 1 , J and Q as a complete set of extensive parameters 
for the mass M(S, J, Q) and define the intensive parameters conjugate to S, Ji and Q. These 
quantities are the temperature, the angular velocities, and the electric potential 

W) Q y fi,= (^) s , c ' * = (w) s y (41) 
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It is a matter of straightforward calculation to show that the intensive quantities calculated 
by eqs. (j39j) - (j4*T]) are the same as those found earlier in this section. Thus, the thermody- 
namic quantities calculated in this section satisfy the first law of thermodynamics, 

k 

dM = TdS + J i + ® d Q. ( 42 ) 

i=l 

B. Stability in the canonical and the grand-canonical ensemble 

The stability of a thermodynamic system with respect to the small variations of the 
thermodynamic coordinates, is usually performed by analyzing the behavior of the entropy 
S(M,Q,3) near equilibrium. The local stability in any ensemble requires that S(M,Q,3) 
be a concave function of its extensive variables or that its Legendre transformation is a 
convex function of the intensive variables. The stability can also be studied by the behavior 
of the energy M(S, Q, J) which should be a convex function of its extensive variable. Thus, 
the local stability can in principle be carried out by finding the determinant of the Hessian 
matrix of M(S,Q,J) with respect to its extensive variables Xi, H^. x . = [d 2 M / dXidXj] 
(33! . In our case the entropy S is a function of the mass, the angular momenta, and the 
charge. The number of thermodynamic variables depends on the ensemble that is used. 
In the canonical ensemble, the charge and the angular momenta are fixed parameters, and 
therefore the positivity of the heat capacity Cj q = T(dS/dT)j^Q is sufficient to ensure local 
stability. The heat capacity Cq,j at constant charge and angular momenta is 

^ (n-1) 

C Q ,j = ^j-[{n - 2)H 2 + l][(n - l)(n - 2)rf ~ 1} + 2q 2 l 2 }{n(n - l)rf ~ 1} - 2q 2 l 2 } (43) 
where T is 

T = 4g 4 / 4 [(3n - 6)H 2 - n + 3)] - 4(n - l)g 2 / 2 r^ 2n ~ 2) [(3n - 6)H 2 - n 2 + 3] 

+n(n - l) 2 (n - 2)r^ 4n " 4) [(n + 2)S 2 - (n + 1)] (44) 

Figure d shows the behavior of the heat capacity as a function of the charge parameter. 
We see that Cqj is positive in various dimensions and goes to zero as q approaches its 
extreme value. Thus, the (n + l)-dimensional asymptotically AdS charged rotating black 
brane is locally stable in the canonical ensemble. In the grand-canonical ensemble, after 
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FIG. 1: Cj 5 q versus q for I = 1, r + = 0.8, n = 6 (bold-line), n = 7 (solid-line), and n = 8 
(dotted-line). 



some algebraic manipulation, we obtain 

1287T 



H 



M 

S,Q,J\ 



n 



[{n - 2)E 2 + l]/ 2 5 6 r 



2^6 T ,3n-4 



n(n - l)r+ (n-1) + 2(2n - 3)Z 2 g 
(n — l)(n — 2)r 



2^2 



2(n-l) 



+ 2/ 2 g 2 



(45) 



As one can see from eq. (|43|l . |H| f Q j| is positive over all phase space. Hence the (n + 1)- 
dimensional asymptotically AdS charged rotating black brane in third order Lovelock gravity 
is locally stable in the grand-canonical ensemble. 



V. CONCLUDING REMARKS 



In this paper, first, we introduced the surface terms for the third order Lovelock gravity 
which make the action well-defined. This is achieved by generalizing the Gibbons-Hawking 
surface term for Einstein gravity or generalizing the surface terms of Gauss-Bonnet gravity. 
We generalized the stress energy momentum tensor of Brown and York in Einstein gravity for 
the third order Lovelock gravity for spacetimes with zero curvature at the boundary. As in 
the case of Einstein gravity, Ig, i& , i& and 1^ of eqs. (JTJ), (JBJ), (JHJ) and (fTT|) are divergent 
when evaluated on the solutions, as is the Hamiltonian and other associated conserved 
quantities. We also introduced a counterterm dependent only on the boundary volume, 
which removed the divergences of the action and conserved quantities of this solution of third 
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order Lovelock gravity. We also found a new class of rotating solutions, whose hypersurfaces 
of constant t and r have zero curvature, in third order Lovelock gravity in the presence of 
cosmological constant and electromagnetic field. These solutions are asymptotically AdS or 
dS for A < or A > respectively. We obtained solutions for special values of «2 and 
given in eq. (J22j) with negative cosmological constant. In the absence of an electromagnetic 
field, these solutions present black branes with one event horizon. The charged solutions 
may be interpreted as black brane solutions with two inner and outer event horizons for 
m > m ext , extreme black holes for m = m ext or naked singularity for m < m ext , where m ext 
is given in eq. (|2fi|l . We found that the Killing vectors are the null generators of the event 
horizon, and therefore the event horizon is a Killing horizon for the stationary solution of 
the third order Lovelock gravity explored in this paper. We computed physical properties 
of the brane such as the temperature, the angular velocity, the entropy, the electric charge 
and potential. Finally, we obtained a Smarr-type formula for the mass of the black brane 
solution as a function of the entropy, the charge and the angular momenta of the black 
brane and investigated the first law of thermodynamics. We found that the conserved and 
thermodynamics quantities satisfy the first law of thermodynamics. We also studied the 
phase behavior of the (n + l)-dimensional charged rotating black branes in third order 
Lovelock gravity and showed that there is no Hawking-Page phase transition in spite of the 
charge and angular momenta of the branes. Indeed, we calculated the heat capacity and 
the determinant of the Hessian matrix of the mass with respect to S, J and Q of the black 
branes and found that they are positive for all the phase space, which means that the brane 
is stable for all the allowed values of the metric parameters discussed in Sec. IIVI This phase 
behavior is commensurate with the fact that there is no Hawking-Page transition for a black 
object whose horizon is diffeomorphic to M? and therefore the system is always in the high 
temperature phase J34|. 
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VI. APPENDIX 



The general solution in n + 1 dimensions is 

f(r) = ^ jl - ( v / 7 + PW + ^)) 1/3 + bj (y 7 + P(r) - fc(r)) V3 J (46) 

where 

, « l « 2 fh-b 2 2 \ 3 . 1 3 1/3 3aA6§ 

and 

w = A+3 ^(^ - („_!)(„ -2^) ) (48) 

The above /(r) reduces to the solution (f2*3*|) when 63 = 63 = !•• At large r the function 
k approaches a constant, and the spacetime is asymptotically de Sitter or anti de Sitter 
depending on the overall sign of 

A eff = ^- |l - (A + v/T + A^) 173 + b 2 2 {-\ + v^TP) 1/3 } 

where we have assumed a > without loss of generality. First, we investigate the conditions 
of the reality of f(r). In order to have real f(r) the expression 7 + k 2 (r) should be positive. 
This occurs for 7 > 0, but for negative 7, this holds if k 2 (r) > The analysis of this case 
proceeds as follows. 1) If A > 0, then k(r) is zero some where and the condition k 2 (r) > \*y\ 
violated, and therefore the function f(r) is complex near the root(s) of k(r) = 0. 2) For 
A < the condition k 2 (r) > \j\ holds provided 

At ^ (\\\ rV~\\ (™- 2 )/ 2 (™- 1 ) / 2 \n/2(n-l) 

A{n — 1) / |A| — a/|7| \ ( q 



(n — 2) \ 30JO3 y \n(n — 1) 

Note that in this case |A| should be larger than a/|t[. Seeking possible black hole solutions, 
we turn to looking for the existence of horizons. The roots of the metric function f(r) are 
located at 

k(r) = k = ± 7 ra (49) 

16 (l + v^l + 4&h 3 J 

which reduces to the value of 1/2 when 7 = as in eq. ()22)1 . Note that ko is real if 
7 3 > -1/(461). Thus > there 

is no black hole solution for 7 3 < —1/ (46?,). For the case of 
uncharged solutions, there exists only one event horizon located at 

m \ V" k - A 

w 



, . a = u , (50) 
2a J 3ati 
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provided a is positive. If a is negative or zero then there will be no roots (corresponding 
to a naked singularity). Note that a = 1/21 2 when 7 = and therefore the above equation 
reduces to eq. (|25|). Again this feature is different from the case of uncharged spherically 
symmetric solutions of third order Lovelock gravity with curved horizon [16j |. for which one 
may have uncharged black holes with two horizons, extreme ones, or a naked singularity. The 
charged solution presents a black hole solution with two inner and outer horizons, provided 
the mass parameter m is greater than m ex t, an extreme black hole for m = m ext and a naked 
singularity otherwise, where m ext is 



Again, one may note that for 7 = (cr = 1/2Z 2 ) eq. (|51|) reduce to m ext given in eq. (|26|). 
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